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ABSTRACT 

Early release science observations of the cluster NGC3603 with the WFC3 on the refurbished HST allow 
us to study its recent star formation history. Our analysis focuses on stars with Ha excess emission, a robust 
indicator of their pre-main sequence (PMS) accreting status. The comparison with theoretical PMS isochrones 
shows that 2/3 of the objects with Ha excess emission have ages from 1 to 10 Myr, with a median value of 3 
Myr, while a surprising 1/3 of them are older than 10 Myr. The study of the spatial distribution of these PMS 
stars allows us to confirm their cluster membership and to statistically separate them from field stars. This result 
establishes unambiguously for the first time that star formation in and around the cluster has been ongoing for 
at least 10-20 Myr, at an apparently increasing rate. 
Subject headings: stars: pre-main sequence - open clusters and associations: individual (NGC 3603) 



1. INTRODUCTION 

The galactic gia nt HII region NGC 3603 located at a dis- 
tance of 7 ± 1 kpc ( iMoffatl 19831) is part of the RCW 57 com- 
plex SE of Eta Car in the Carina arm. The young, bright, 
compact stellar cluster (HD97950 or NGC 3603 YC) lies at the 
core of this region and has long been the center of attention for 
the relatively numerous population of massive stars at its cen- 
ter. The collection of 3 WNL, 6 03 and numerous late O type 
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stars together with a bolometric luminosity of 100 times that 
of the Orion cluster and 0. 1 times that of NGC 2070 in the 30 
Dor complex in the Large Magellanic Cloud (LMC) and a to- 
tal mass in excess of 10 4 M Q (Harayama et al. 2008, hereafter 
HA08, and references therein) places it squarely in the cate- 
gory of a super star cluster of the type more commonly seen in 
star burst regions of young galaxies. Thus, it represents a re- 
solved possible pr ototype of the building blocks of forming or 
merging galaxies (lLarsenll2009l) and an excellent opportunity 
to better understand their complex star formation history and 
the physical processes involved (see, for example, P eters et al. 
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In this context, NGC 3603YC (from now on referred to 
as simply NGC 3603) certainly contains many stars younger 
than 3 Myr but it is so far unclear if it also contains older 
stars as might be sug gested by the prese nce in the field of the 
evolved star Sher25 (Mel enaet al.|2008l) , located ~ 20" north 
to the cluster center. This exciting possibility has yet to be es- 
tablished unambiguously because the latter's membership in 
the cluster is still quite uncertain. The determination of the 
cluster's initial mass function (IMF), of course, also depends 
critically on its formation history as imprinted in the age dis- 
tribution that can only be determined accurately once we dis- 
entangle the various populations from each other. Results on 
this c lusterOs IMF (HA08: IStoite et alJl2006tlSung & Besselll 
2004) may be affected by imperfect correction for this effect. 
In order to investigate the star formation history in 
NGC 3603, we have used the Wide Field Camera 3 (WFC3) 
on Hubble Space Telescope (HST) in the UVIS mode to im- 
age NGC 3603 with broad band and narrow band filters as 
part of the early release science (ERS) program. In this pa- 
per, we describe the first results of this investigation. The ob- 
servations and analysis procedure are presented in Section [2] 
In Section [3] we show the resulting color magnitude diagram 
(CMD) and the estimation of field contamination and extinc- 
tion. In Section 2] we describe the identification of pre-Main 
Sequence (PMS) stars, while their spatial distribution and a 
discussion of the preliminary physical implications of the re- 
sults can be found in Section[3]and[6] respectively. 
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FIG. 1 . — WFC3 I'.l X 2',1 mosaic images in the F555W (left panel) and F656N (Ha line; right panel) filters of NGC 3603 star forming region as obtained 
through PyRAF/MULTIDRIZZLE package. North is 30° to the right of the vertical. East is to the left of North. 



2. OBSERVATIONS AND DATA REDUCTION 

The photometric data used in this work consist of a series 
of deep multi-band images acquired with the new WFC3 on 
board of the HST. The WFC3 consists of two detectors, one 
optimized for observations in the wavelength range ~ 200 to 
~ 1000 nm (UVIS channel) and the other between ^0.9 and 
~ 1.7 fim (IR channel). The UVIS detector consists of two 
2Kx4K CCDs covering a field of view (FoV) of 1 62"x 1 62" at 
a plate scale of 0."04/px. The IR detector is a single lKxlK 
HgCdTe CCD offering a total FoV of 123"x 136" at a pixel 
resolution of 0"13. A more detailed description o f the WFC3 
and its current performances can be found in Wong et al. 
(2010). 

The data used in this work are part of the ERS observa- 
tionQ obtained by the WFC3 Scientific Oversight Commit- 
tee for the study of star forming regions in nearby galaxies 
(Program ID number 1 1360). NGC 3603 was observed using 
both the UVIS and IR channels. In this paper we use the im- 
ages taken through the broad-band F555W and F814W and 
narrow -band F656N filters for a total exposure time of 1000s, 
1550s and 990s, respectively. The IR data-set will be pre- 
sented in a separate paper (Spezzi et al. 2010, in preparation). 

Three images with approximately the same exposure times 
were taken with a few pixel dithering in order to allow for the 
removal of cosmic rays, hot pixels and other detector blem- 
ishes. All the observations were performed so that the core 
of NGC 3603 is roughly at the center of the camera's FoV. In 
Figure [1] we show a mosaic of the images in the F555W (left 
panel) and F656N (right panel) filters as obtained with the 
PyRAF/MULTIDRIZZLE package. 

The photometric analysis of the entire data-set was per- 
formed on the flat-fielded (FLT) images by adopting the fol- 
lowing strategy. The images, corrected for bias and Hatfield, 
need a further field-dependent correction factor to achieve 

23 A complete description of the program, targets and the observations can 
be found at the web page http://archive.stsci.edu/prepds/wfc3ers/ 



uniformity in the measured counts of an object across the 
field. Applying the correction simply involves multiplying the 
FLT images by the pixel area map images. A large number of 
isolated, well exposed stars were selected in every image over 
the entire FoV in order to properly model the point spread 
function (PSF) with the DAOPHOTII/PSF routine dStetsonl 
1987). We used a Gaussian analytic function and a second or- 
der look-up table was necessary in order to properly account 
for the spatial variation in the images. 

A first list of stars was generated by searching for ob- 
jects above the 3a detection limit in each individual im- 
age and a preliminary PSF fitting run was performed using 
DAOPHOTII/ALLSTAR. We then used DAOMATCH and 
DAOMASTER to match all stars in each chip, regardless of 
the filter, in order to get an accurate coordinate transformation 
between the frames. A master star list was created using stars 
detected in the F814W band (the deepest in the UVIS data- 
set) with the requirement that a star had to be detected in at 
least two of the three images in this filter. 

We used the sharpness (sh) and chi square parameters given 
by ALLSTAR in order to re move spurious de t ections. It has 
alread y been shown (see e.g. ICool et al.lll996t lAscenso et al.l 
120071) that these parameters are good tracers of the photo- 
metric quality. By using a sample of real stars we found 
the range -0.15 < sh < 0.15 to be safe enough to eliminate 
most spurious objects. The final catalogue was then obtained 
by rejecting any residual spurious source (mostly associated 
with bright emission peaks in the HII region not due to point 
sources) through visual inspection of the drizzled imag es. The 
maste r list was then used as input for ALLFRAME ( Stetson 
1994), which simultaneously determines the brightness of the 
stars in all frames while enforcing one set of centroids and 
one transformation between all images. All the magnitudes 
for each star were normalized to a reference frame and aver- 
aged together, and the photometric error was derived as the 
standard deviation of the repeated measurements. The final 
catalogue of the UVIS F555W, F656N and F814W bands 
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FIG. 2. — Color magnitude diagram of NGC 3603 as obtained from WFC3 
observations. The ZAMS from Marigo et al. 1 2008) at the solar metallicity 
for an assumed distance modulus (tn-M)o = 13.9 and average extinction Ay = 
4.5 is over-plotted (dashed line). The average photometric magnitudes and 
color errors are shown (crosses). The reddening vector (Av=l) is also shown. 

contains around 10,000 stars. 

The photometric c alibration was performed follow- 
ing |KajjraTer]n](|2(J09|). A sample of bright isolated stars was 
used to transform the instrumental magnitudes to a fixed aper- 
ture of 0."4. The magnitudes were then transformed into the 
VEGAMAG system by adopti ng the synthetic zer o points for 
the UVIS bands (see Table 5 of lKalirai et aT]|2009l) . Hereafter 
we will refer to the calibrated magnitudes as V, I , and Ha to 
indicate nifmwi mpzuw and nif^6N, respectively. 

The WFC3/UVIS channel is affected by geometric dis- 
tortion and a correction is necessary in order to prop- 
erly derive the absolute positions of individual stars in 
ea ch catalogue. We used t he dis tortion coefficients derived 
bv lKozhuri na-Platais et alj ( 120091) to obtain relative star posi- 
tions that are corrected for distortion. We then used the stars 
in common between our UVIS and the 2MASS catalogues to 
derive an astrometric solution and obtain the absolute RA and 
DEC positions of our stars. We find a systematic residual of 
~ 0"3 with respect to the 2MASS coordinates. 

3. THE UVIS COLOR MAGNITUDE DIAGRAM 

The CMD as obtained with our photometric reduction pro- 
cedure is shown in Figure[2] Typical photometric errors (mag- 
nitudes and colors) are indicated by black crosses. This is 
the deepest and most accurate optical diagram derived so far 
for this cluster, demonstrating the unprecedented photometric 
capabilities and the resolving power of the new WFC3. We 
show in the CMD t he pos ition of the zero age MS (ZAMS) 
from Marigo et al. (2008) for solar metallicity (solid line), 
having adopted a distance modulus (m-M)o =13.9 (see HA08 
and references therein), the average value of the extinction 
Ay = 4.5 as reported in the literature (see ISung & Besselll 
12004 hereafter SB04), and having assumed the extinction law 
of lCardelli et al] d!989h . Note that the value of A v = 4.5 is as- 
sumed here only for illustration purposes and for comparison 



with previous studies, as we will show later that the actual ex- 
tinction value is larger in the area that we studied. Nonethe- 
less, the position of the ZAMS obtained in this way helps us to 
identify a population of candidate MS stars that extend from 
the saturation li mit a t V ~ 1 7 down to V ~ 26. As we will 
show in Section [3~T1 and l3~2l most of these objects are likely 
foreground field stars but some of them represent a bona-fide 
low mass MS population belonging to the cluster. 

A discontinuity in the stellar color distribution in this CMD 
separates a population of objects along the ZAMS from one 
clearly grouped at redder colors and consistent with the young 
population of PMS stars already detected in NGC 3603 (see 
e.g. HA08; SB04). This confirms that NGC 3603 is an active 
star forming region. In order to learn more about the proper- 
ties of this recent star formation episode, it is useful to com- 
pare our CMD with PMS isochrones. However, this requires 
detailed knowledge of possible sources of error such as the 
contamination from field stars and the presence of differential 
reddening. We address these issue here below. 

3.1. Field star contamination and extinction 

As already discus sed by HA08 

and |Nurnberger & Petr- Gotzensl (120021) . the region of 
the CMD fitted by the ZAMS in our CMD is contaminated by 
field stars with a luminosity distribution roughly in agreement 
with Galactic models. In order to quantify the degree of 
contamination, we generated a simulated catalogue of field 
stars using the Galactic model of Robin et al. (2003). The cat- 
alogue covers a projected area of 162" x 162", corresponding 
to the FoV of the WFC3 (see Section EJ), with a diffuse 
extinction of 0.7mag/kpc in a distance interval of 7kpc (i.e. 
the cluster distance, see Section Q]). Since the magnitudes 
of the synthetic stars are given in the Johnson-Cousins (JC) 
photometric system, we used model atmospheres from the 
ATLAS9 library of Kurucz ( 1993) to calculate the magnitude 
difference between the JC and WFC3 photometric systems 
as a function of the effective temperature of the stars. Com- 
parison of the synthetic catalogue with our CMD reveals that 
about 85% of the stars in the region around the ZAMS in our 
photometry are potentially field stars, while it decreases to a 
few percentage towards the region populated by PMS stars. 

Recently, Rochau et alJ (2010) presented a proper motion 
study of NGC 3603 based on HST/WFPC2 (Wide Field Plan- 
etary Camera 2) observations obtained ten years apart, respec- 
tively in 1997 and 2007. On the left panel of their Figure 2 
these authors show the position of field stars (open circles) on 
the CMD obtained with the Planetary Camera (PC) data. By 
comparing this diagram with the same obtained when con- 
sidering only bona-fide cluster members (i.e. those with a 
similar proper motion; central panel on the same figure), we 
quantify the field contamination along the MS in the PC data 
to be of order 50%. Although not explicitly stated in their 
paper, also Bra ndl et al] ([1999) reach a similar conclusion, as 
about 50 % of the stars on their original lower MS are still 
present after statistical subtraction of a comparison field (see 
their Figur e 2). Con sidering that both the PC dataset used 
by Rochau et al. (2010) and the observations of Brandl et al. 
(1999) sample the core region of the cluster, where the den- 
sity of cluster stars is highest, we conclude that the 85 % value 
of field star contamination that we find in th e external re- 
gions from the models of iRobin etall d2003l) is reasonable 
and might actually be an overestimate of the true value. We 
therefore assume it as an upper limit to the contamination 
level in this field. We will return in Section [3~2l to the work of 
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iRochau et all (120101) and iBrandl et ail (1 1999b to discuss their 
implications for the study of the cluster's stellar population, 
but we first need to address the issue of differential extinction 
in this field. 

An efficient method to quantify the amount of differential 
extinction in a star cluster is to use the position of a star in the 
observed CMD to calculate its distance from a fiducial line 
(e.g. the ZAMS itself) along the direction of the reddening 
vector. This distance would be the resultant of two compo- 
nents, namely E(V-I) on the abscissa and Ay on the ordinate, 
and would give us an esti mate of the extin ction toward the star 
itself (see an example in Piott o et afll 19991) . 

This method works under the assumption that the observed 
stars are at the same distance, i.e. belong to the same system. 
As discussed above, the field star contamination in the range 
of magnitudes covered by our data is high. The described 
method can be applied in a reliable way only to bright ob- 
jects (V< 17), since these are very likely cluster members and 
field star contamination is minimized at these magnitudes. Al- 
though our data cannot be used for this purpose, since all stars 
brighter than V ~ 17.5 are saturated, using shorter exposures 
SB04 were able to perform a study of differential reddening 
in NGC 3603 taking advantage of multi-band HST photome- 
try of the bright massive stars (i.e. the same objects that are 
saturated in our images). These authors were able to map 
the variation of E(B- V) as a function of the distance from the 
cluster centre (see their Figure 5b). They found that the value 
Ay ~ 4.5 is representative of the very centre of the OB stars 
association, while they noticed an increase toward the external 
regions. Following the work of SB04 and adopting Ry = 3.55 
as they suggest, we estimate that the mean value of Ay in the 
area sampled by our observations (from ~ 10" to ~ 70") is 
Ay = 5.5. Therefore, in the rest of this paper we will adopt 
this value to correct our magnitudes for extinction. The final 
CMD corrected in this way is shown in Figure[3] 

3.2. The PMS population in NGC 3603 

Once extinction is taken into account, the CMD shown in 
Figure |3] can be used to determine stellar ages through PMS 
isochrone fitting. PMS isochrones with ages of 1, 2, 3, 10, 
20 and 30 Myr from Sies s et al.l (|2000) are shown in the fig- 
ure from right to left. We used the same value of t he distance 
modul us adopted for the ZAMS fit. Note that the ISiess et all 
(2000) models are not available for the WFC3 phot ometric 
system. As above, we used the ATLAS9 library of iKuruczl 
( 1993) to calculate the magnitude differences between the JC 
to WFC3 photometric systems. After this correction, the mod- 
els indicate that the lowest mass that we reach for PMS stars 
is0.3M Q . 

In the magnitude range V < 20, where the photometric un- 
certainty on the V— I color (~ 0.05) is smaller than the typical 
isochrone separation, the CMD suggests for our PMS stars an 
age in the range from 1 to 10 Myr, with an average age of 
3 Myr. Assuming the PMS isochrones are correct, this can 
already be considered as tentative evidence of an age spread 
in the stellar population of NGC 3603. 

Recently, HA08 studied the IMF of NGC 3603 using near- 
infrared (IR) imaging from ground-based adaptive optics 
photometry of the cluster center obtained with the NAOS- 
CONICA (NACO) camera and the wider field Infrared Spec- 
trometer and Array Camera (ISAAC) at the Very Large Tele- 
scope (VLT). Their JHKL 1 bands photometry reaches the 
magnitude limit J ~ 20.5 (i.e. ~ 0.4M Q for cluster stars) in 
an area of ~ 110" radius from the cluster center. By com- 
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FIG. 3 .— De-reddened CMD. The ZAMS from (Marigo et al. 2008, dashed 
line) is used to identify the location of MS stars. The population of PMS stars 
is fitted with 1, 2, 3, 5, 10, 20 and 30 Myr PMS isochrones (from right to the 
left) from Siess et al. (2000) (solid lines). The position of the 1 and 10 Myr 
isochrones is indicated. According to the models a mass of M ~ 0.3Mq is 
reached at magnitude V= 23. 



paring their CMD to lBaraffe et all (1998) PMS evolutionary 
models, they identify a population of PMS stars with ages of 
0.5- 1.0 Myr. Moreover, by adopting a set of MS isochrones 
from Lejeune & Schaerer (2001), they provide tentative hints 
of the presence of an evolved MS population of 2. 0-2. 5 Myr. 
Note that this age estimate is based on the comparison of the 
MS isochrones with the position of three massive evolved 
O stars in their IR CMDs. Coupling this piece of evidence 
with the presence of the evolved post-red supergiant star 
Sher25 dMoffat 19831) in the cluster field, HA08 hypothe- 
size a possible age spread in the cluster population suggest- 
ing the presence of two distinct bursts in the star formation 
history, separated by ~ 10 M yr. It is interesting to note that 
recentlv livlelena et al.l J200 8) have placed Sher 25 at the same 
distance as NGC 3603, suggesting a common origin. 

The same ISAAC observations analyzed by HA08, were 
previously used by Stolte et al. (2004) and Brandl et al. 
( 1999) to study the low-mass stars population in NGC 3603. 
While by inspecting the CMDs shown in their Figure 4 and 
3c, respectively, signatures of the presence of a low MS pop- 
ulation in the CMDs can not be ruled out, both these papers 
agree on dating a PMS population in the range of 0.5- 1 Myr. 

SB04 published what was by then the deepest optical 
CMDs based on a combination of UBVI and Ha photome- 
try from the Siding Spring Observatory (SSO) and archival 
HST/WFPC2 observations. The SSO ground-based CMDs 
sample the brightest massive population. The WFPC2 CMD 
samples stars down to V ~ 21. In the very central region of 
the cluster (r < 12"), the CMD shows the presence of a well 
defined population of massive MS stars together with a pop- 
ulation of low mass PMS objects. Apparently only few MS 
stars are detected in the cluster center at magnitudes fainter 
than V ~ 19.5, while in the external regions (18"- 120") stars 
with masses down to ~ 1M© are detected both in the MS and 
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PMS regions. 

By using archival WFPC2 observations in the F656N band 
combined with the WFPC2 broad-band photometry, SB04 
were able to identify objects with excess Ha emission, which 
is a signature of the PMS phase. While the majority of 
these objects on their CMD occupies a region consistent with 
very young PMS isochrones (1 Myr), some stars fall near the 
ZAMS (see their Figure 7). SB04 speculate on a possible 
spread in age of the cluster stars, but they warn that it is dif- 
ficult to reach a firm conclusion because of the decreasing 
completeness and photometric accuracy of their photometry 
at fainter magnitudes (V > 19). 

A careful examination and comparison of the recent stud- 
ies on NGC 3603 mentioned above shows that they do not 
exclude the possible presence of multiple star formation 
episodes in the recent cluster history and even of a popula- 
tion of low mass MS stars. However, this scenario necessarily 
lacks a clear observational confirmation, since field contami- 
nation remains a crucial point in the study of stellar popula- 
tions in NGC 3603. An efficient way to overcome these ob- 
stacles, when observations exist that are separated by a suffi- 
ciently large temporal baseline (> lOyr), is the use of proper 
motions to separate cluster stars from foreground and back- 
ground objects. 

As mentioned in Section lXTl the proper motions s tudy for 
the core region of NGC 3603 bv lRochau eTaTI (120 1 Ot) showed 
signatures of at least two star formation epochs in the cluster, 
1-2 Myr and 4-5 Myr old, respectively. Moreover, the CMD 
as cleaned-up from the field stars, reveals the presence of a 
low mass MS. The authors conclude that the latter is likely 
a population of objects not belonging to the cluster, that the 
proper motion technique failed to identify. We will offer later 
in our paper new observational evidence supporting the hy- 
pothesis that an old stellar population belonging to the cluster 
is present in the region of t he CMD occupied b y the candidate 
low MS stars identified bv lRochau et all (120101) . 

As an alternative and independent method to investigate the 
star formation history in NGC 3603, we have decided to take 
full advantage of our new deep F656N band exposures to 
search for objects with excess Ha emission, since this feature 
is a good indicator of the PMS stage and therefore of recent 
star formation. By looking at the spatial distribution and age 
of all the objects with Ha excess emission, we will be able to 
better understand their cluster membership. 

4. Ha EMISSION STARS 

The presence of a strong Ha emission line (EW# a >10 A) 
in young stellar objects is normally interpreted as a sig- 
nature of the mass accretion process onto the surface of 
the object that requires the presence of an inner disk (see 
iFeigelson & Montmerld [19991 IWhite & BasrHHOOl and ref- 
erence therein). 

The traditional approach to search photometrically for Ha 
emitters is based on the use of the R-band magnitude as a 
measure of the level of the photospheric continuum near the 
Ha line, so that stars with strong Ha emis sion will have a 
largei? — Ha color. However, as discussed in De Marc hTet al.l 
(2010a), since the R band is over an order of magnitude wider 
than the Ha filter, the R-Ha color does not provide an accu- 
rate measurement of the stellar continuum level inside the Ha 
band. Thus, while helpful to identify PMS stars, the R - Ha 
color does not provide an absolute measure of the Ha lumi- 
nosity nor of the Ha equivalent width. This additional in- 
formation can be derived using measurements in the neigh- 
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FIG. 4. — Selection of Ha emission excess stars in a color-color diagram. 
The dashed line represents the median V — Ha color representative of stars 
with no Ha excess. Solid triangles and circles show the position on the dia- 
gram of the counterpart of the 67 objects with Ha excess emission by SB04 
accepted and rejected through our selection criteria, respectively. 



boring V and I bands, as recently shown by De Marc hTet alJ 
(I2010al) . This method allows us to reliably identify PMS ob- 
jects actively undergoing mass accretion regardless of their 
age. Briefly, the method combines V and I broad-band pho- 
tometry with narrow-band Ha imaging to identify all stars 
with excess Ha emission an d to measure their Ha lu minosity 
and mass accretion rate (see lDe Marchi et aTJl2010al for more 
details). 

We followed this approach to select bona-fide PMS stars in 
the field of NGC 3603. Figure g] shows the V-Ha vs. V-I 
diagram for the stars in our catalogue. We use the median 
V-Ha dereddened color of stars with small (< 0.05 mag) 
photometric uncertainties in each of the three V, I and Ha 
bands, as a function of V —I, to define the reference template 
with respect to which the excess Ha emission is identified 
(dashed line in Figure |4). We selected a first sample of stars 
with excess Ha emission by considering all those with a V - 
Ha color at least 5 a above that of the reference line, where a 
here is the uncertainty on the V - Ha color of the star. Then 
we calculated the equivalent width of the Ha emission line 
(EWgq) from the measu red color excess using Equation 4 of 
iDe Marchi et al.l d2010 a). We finally considered as bona-fide 
PMS stars those objects with EW Hce >10 A dWhite & Basril 
2003) and V-I > 0; this allows us to clean our sample from 
possible contaminants, such as older stars with chromos pheric 
activity and Ae/Be stars, respectively (see lScholz et alj 2007). 

With this approach, we selected a first sample of ~ 800 ob- 
jects with Ha excess emission. Through a visual inspection 
of the images, we noticed that some of these objects, although 
well detected both in the V and I bands, are located along fil- 
aments of gas and dust clearly visible in the Ha image (see 
right panel in Figure [TJ. It is crucial to consider that, if the 
centroid of a star falls on top of a filament that is only partially 
included in the annulus that our photometry routines use for 
background subtraction (from 4 to 7 pixel radius), the back- 
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FIG. 5. — Four regions (each about 9" x 5" in size) of the drizzled Ho image as obtained by applying an unsharp-masking algorithm to highlight and sharpen 
the details of the dusty and cloudy structures. All candidate Ho excess emitters are marked with two concentric annuli (4 and 7 pixels radius, respectively). The 
objects excluded from our bona-fide sample are marked with a cross. 



ground will be underestimated and the derived Ha magnitude 
of the star will be over estimated. All sources with Ha excess 
emission have been visually inspected and all those falling on 
top of one of such filaments conservatively excluded. 

As an example, we show in Figure [5] four regions (each 
about 9" x 5" in size) of the drizzled Halpha image (with the 
highest spatial resolution and free of cosmic rays) as obtained 
by applying an unsharp-masking algorithm to highlight and 
sharpen the details of the gas filaments. All stars with ap- 
parent Ha excess emission are marked on the Figure with 
two concentric annuli corresponding to the area in which the 
background has been estimated (4 and 7 pixel radius, respec- 
tively). We have marked as suspicious and excluded from our 
bona-fide sample all those with significant and non-uniform 
filament contamination inside the photometric aperture (see 
objects marked with a cross in Figure |5}. In this way, we re- 
duced our sample to 412 objects with 5a Ha excess emission 
that we consider bona-fide PMS stars. 

In light of the considerable amount of differential redden- 
ing present in our field, one issue that needs to be taken into 
account is the impact that extinction will have on our selec- 
tion of objects with Ha excess emission. Uncertainties on the 
extinction will move the objects in Figure |4] mostly along the 
V—I axis, thereby changing the reference value of the V—Ha 
color. In other words, variable extinction would change the 
reference relation (dashed line) to redder or bluer colors for 
individual stars. In order to estimate the impact of differential 
reddening, we have simulated the uncertainty introduced by a 
variation of ±0.2 mag in E(V— I), According to SB04, such a 
value is the characteristic reddening variation seen in the area 
covered by our observations (see their Figure 5). Furthermore, 
it corresponds to AAv = ±0.5 mag, which fully covers this 
spread between the Av = 5.5 value that we find for the redden- 
ing and the canonical figure for NGC 3603 (Av = 4.5). If the 
true E(V -I) value of a star were underestimated, we would 
also underestimate its V — Ha excess, thereby in practice im- 
posing a more stringent limit on the excess value itself. By 
underestimating E(V -I) by 0.2 mag , our 5a detection limit 
would correspond to about 6 a. Conversely, by overestimating 
the true E(V -I) by 0.2 mag, about 17 % of the stars detected 



at the 5 a level will drop below that, but all of them would still 
be above the A a level. Therefore, while there might be small 
uncertainties as to the value of the Ha excess for a specific 
object, the adoption of an average extinction value across the 
field does not significantly affect our selection of bona-fide 
PMS stars. 

In Figure[6]we indicate with black dots the positions in the 
average dereddened CMD of all objects with excess Ha emis- 
sion. The comparis on with the theoretical 1, 3 and 10 Myr 
PMS isochrones of ISiess et"ai . (2000) (dotted-dashed, solid 
and short-dashed line, respectively) suggests a typical age of 
~ 3Myr for 2/3 of the objects with Ha excess emission, but 
a surprising 1 /3 of them are lying at or near the ZAMS (long 
dashed line), thus suggesting a considerably higher age, of or- 
der 20-30 Myr. A 10 Myr PMS isochrone seems to efficiently 
divide the stars with Ha excess emission in two populations, 
one lying near the ZAMS and the other in the PMS region. 

Even though the fraction of PMS stars showing accr e tion is 
known to drop rapidly with age, ISicilia-Aguilar at ail (120 lOt) 
recently found evidence of stars older than 10 Myr still under- 
going mass accretion (see also Nauven et al. 2009, and refer- 
ences therein). In the LMC, De March i et al.1 (I2010al) found 
about 130 actively accreting PMS stars with a median age of 
14 Myr indicating that ~ 10 Myr old stars can show Ha excess 
emission due to the accretion process. 

However, as already discussed in Section |3~T1 the region of 
the CMD fitted by the ZAMS track is strongly contaminated 
by field stars. It is then crucial to understand whether the ob- 
jects with Ha excess emission that we see in NGC 3603 are 
indeed cluster members belonging to a previous generation 
or are just very young field stars. We will show in Section [5] 
the results of a statistical analysis of stellar membership pro- 
viding convincing evidence that these older PMS stars have a 
different distribution from that of field objects. 

It is known that one property of the accretion process that 
characterizes PMS stars is temporal variability, which can be 
traced to changes both in the continuum and more importantly 
in the emission lines (e.g. Herbst 1986; Harti gan et al. 1991; 
Nguyen et al. 2009). As mentioned in SectionEU SB04 pub- 
lished a list of sources with Ha emission in NGC 3603 de- 
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FIG. 6. — The location of stars with Ha excess (black points) are shown on 
the CMD together with a ZAMS and the 1, 3 and 10 Myr PMS isochrones 
adopted in Figuref5](long dashed, dotted-dashed, solid and short-dashed line, 
respectively). As expected, most of the selected objects are young PMS stars 
(age < 10 Myr; YPMS box). Surprisingly, a number of sources are located on 
or near the ZAMS isochrone, i.e. in a position of the CMD compatible with 
a older PMS population (age > 10 Myr; OPMS box). The selection boxes of 
the two populations are shown (black lines). 

rived from SSO and WFPC2 observations. We searched for 
the counterparts of these objects in our observations by cross- 
correlating our photometric catalogue with that of SB04. The 
stars detected by SB04 on the SSO data are bright and are all 
saturated in our images. On the other hand, the WFPC2 cat- 
alogue of SB04 contains 96 sources selected on the basis of 
their Ha index, of which 67 have a counterpart in our cata- 
logue. As for the 29 missing objects, ten fall in the inner 10" 
of the cluster center, where saturation in our images makes 
star detection very difficult. Two of the remaining 19 sources 
fall outside the UVIS FoV, while the others are very close to 
the saturation limit in our catalogue or fall into the gaps of the 
UVIS mosaic. 

All these 67 sources had Ha excess emission at the time of 
the SB04 observations, but only 35 of them have Ha excess 
emission at the > 5 a level when our observations were taken. 
This number decreases to 23 if we consider only stars with 
EW{f a > 10 A. When lowering the threshold to 3 a, the num- 
ber of stars with Ha excess in common with SB04 grows to 
41, but only 25 of them have EW// Q >10 A. This means that 
most (~ 65 %) of the sources showing Ha excess in the study 
of SB04 (March 1999) do not show it at the epoch of our ob- 
servations (Aug 2009). This is in line with the current under- 
standing of the accretion mechanism whereby the in-falling 
material from the circumstellar disk is subje ct to bursts, cor- 
respo nding to peaks in the Ha emission (e.g. Fernandez et al. 
1995). 

SB04 looked for matches between bright X-ray sources and 
objects with Ha excess emission by using archive Chandra 
X-Ray Observatory images. Only ~ 1/3 of the 96 Ha excess 
emission sources in the WFPC2 data set have a X-ray emis- 
sion. We found that of the 23 stars with Ha excess emission in 



common between our catalogue and that of SB04, a total of 8 
have X-ray emission, according to SB04. While it is expected 
PMS stars actively undergoing mass accretion to be detected 
in X-rays, it is also true that the physical processes related to 
X-ray flaring and those pertaining to mass accretion are not 
the same (see e.g. Feigelson 2005). As shown by SB04, we 
should not expect a complete match between the sources that 
are bright in X rays and those showing Ha excess at any given 
time. And, obviously, even less so when the observations are 
not simultaneous, as in the case of the Chandra and HST data 
used by there authors, because of the considerable variability 
to which the accretion process is subject (see above). 

5. SPATIAL DISTRIBUTION OF PMS STARS 

The presence of stars with Ha excess emission overlap- 
ping the ZAMS in the CMD of Figure [6] offers new support 
to the hypothesis of a spread in the age of the PMS stars in 
NGC 3603 (see HA08 and references therein), since this is 
the region where PMS stars older than 10 Myr are expected 
to be located. One obvious issue to address is whether these 
objects are cluster members, as we know that there is a po- 
tentially sig nifica nt contribution from field stars in this region 
(see Section l3~Tl i. In order to investigate cluster membership 
for stars showing Ha emission, we can look at their spatial 
distribution compared to that of field stars. 

Based on the distribution of Ha excess objects in the CMD, 
we define two regions containing bona-fide PMS stars (i.e. 
those with Ha excess of different ages). The 10 Myr PMS 
isochrone is used as a guide to define a rough separation 
between the young PMS population (< 10 Myr; hereafter 
YPMS) and the old PMS population (> 10 Myr; hereafter 
OPMS). We finally define as field population all stars lying 
in the OPMS box and not showing Ha excess. The corre- 
sponding selection areas are shown as boxes in Figure [6] It is 
important to underline here that, having defined as bona-fide 
PMS stars only the objects with excess Ha emission, we are 
in practice setting a lower limit to the actual number of stars 
in the PMS phase in NGC 3603, since some of them can be in 
the PMS stage without showing any Ha excess because of the 
variability of their Ha flux (see Section HJ. This is certainly 
the case for the largest majority of stars that lie to the right of 
the 10 Myr isochrone in Figure 3. On the other hand, since we 
are interested in the radial distribution of PMS stars, select- 
ing only objects with Ha excess emission does not affect the 
significance of our statistical analysis. 

The cumulative radial distribution of the four groups de- 
fined above (PMS, YPMS, OPMS and Field) with respect 
to the cluster center is shown in the left panel of Figure [7] 
whereas in the right panel the positions of YPMS stars 
(crosses) and of OPMS stars (filled circles) are shown on the 
F656N band image. To calculate the radial distribution of 
these objects we adopted the RA(J2000) = ll h 15 m 7. 8 26 and 
DEC(J2000) = -61°15'37."9 as coordinates for the cluster 
center, following SB04. We exclude from our analysis the 
innermost 5" radius where there is a high concentration of 
bright O-B type stars that are saturated in our long exposures, 
and the high crowding level makes object detection difficult. 
The radius r = 70" defines the largest circle inscribed in the 
FoV. 

The graph on the left panel of Figure [7] clearly shows that 
PMS stars (solid line) are more centrally concentrated than 
field stars (dot-dashed line). Furthermore, the radial distri- 
butions of YPMS and OPMS stars (dotted and dashed lines, 
respectively) are clearly different from the field population, 
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FIG. 7. — Left panel-Cumulative radial distribution of the four star groups as indicated in the legend. Right panel- The location of YPMS (crosses) and OPMS 
(solid circles) over-imposed on the F656N image already shown on Figure[T] 



supporting the idea that these stars belong to NGC 3603. We 
used a Kolmogorov-Smirnov (KS) test to check the statisti- 
cal significance of the differences in the observed distribu- 
tions. The test yields more than 3<r confidence level that field 
stars have a different radial distribution from that of the PMS, 
YPMS and OPMS groups. 

Interestingly, YPMS stars appear to be more centrally con- 
centrated than OPMS objects, contrary to what one would 
expect in a triggered star formation scenario. Note that the 
KS test, as used here, indicates the probability that the differ- 
ent groups are drawn from the same population on the basis 
of their radial distribution with respect to a common center. 
In principle, however, the assumption of a common center of 
gravity for the PMS population as a whole could be incorrect, 
since the center of gravity of the OPMS stars might as well 
differ from that of the younger YPMS population. Nonethe- 
less, the analysis presented here clearly shows how powerful 
the use of the Ha excess information can be to identify bona- 
fide PMS stars and to properly separate them from the field 
star population. 

6. DISCUSSION AND CONCLUSIONS: THE STAR FORMATION 
HISTORY OF NGC 3603 

The literature on the study of the stellar population of 
NGC 3603 offers a wide debate on the age of the st ars in 
this cl uster and on their formation history. While Stolte et al. 
(2004) and SB04 suggest 1 Myr as a common age for 
massive MS stars and PMS objects, possible evidence of 
larger age spreads comes from the analysis of the popula- 
tion of massive stars. HA08 derive an age of 2 -2.5 Myr 
for MS stars through isochrone fitting of three massive stars 
in the center of the cluster classified as WN6h+abs objects 
by lCrowther & Dessa rt ( 1998), while they date the PMS stars 
in the range 0.5- 1.0 Myr. Hendry et al. (2008) also notice 
that the massive-star population in NGC 3603 appears to be 
predominantly coeval (with an age of 1 - 2 Myr), but that 
Sher25, and one O-type supergiant, are likely to be slightly 
older (-4-5 Myr; Melenaet al. 2008; Crowtheret all 12006). 



Sher25 (Sher 1965; lMoffatl[T983l) is clearly visible in Fig- 
ure Q] and has been recently classified as a blue supergiant 
(BSG) surrounded b y an asymme t ric, h ourglass-shaped cir- 
cumstellar nebula by Hendry et al. (2008). 

HA08 suggests the possibility that, if Sher 25 belongs to 
the cluster, two distinct star formation episodes separated by 
~ 10 Myr should be present in NGC 3603. Sher 25 is not the 
only B SG in the NGC 3603 region. Spectroscopy by Moffat 
( 1983) revealed two other B SGs in the vicin ity of the cluster 
core (see also Figure 1 from|Brandner et al. 1997). The latter 
authors conclude that the simultaneous presence of BSGs and 
stars of MK type 03 V requires at least two distinct episodes 
of star formation in NGC 3603 separated by ~ 10 Myr. 

Our study shows that the OPMS population belongs to 
NGC 3603 and not to the field, thus establishing unambigu- 
ously for the first time that the star formation of this cluster 
is not characterized by a 1-2 Myr single burst. These objects 
might indeed represent the low mass population of a star for- 
mation episode that occurred more than 10 Myr ago, which 
would naturally explain the presence in the same region of 
evolved massive stars (age > 10 Myr) like Sher 25. 

In order to better understand how star formation has pro- 
ceeded recently in NGC 3603, we can look at the age distri- 
bution of PMS stars in the CMD. Shown in the left panel of 
Figure|8]are the PMS isochrones of Siess et all (120001) for ages 
of 1, 2, 4, 8, 16, 32 and 64 Myr. Note that the photometric un- 
certainty of our data (see large crosses on the figure) confirms 
that it is possible to assign relative ages to these stars with 
an accuracy of better than a factor of 2. Solid thick points 
and asterisks in the figure represent all bona-fide PMS objects 
(i.e. those with Ha excess emission) having ages in the range 
1-64 Myr and masses between 0.9 and 4M Q . The corre- 
sponding age distribution is shown by the solid histogram in 
the right panel of the figure. In the same panel we also show, 
as a dashed line, the histogram of the age distribution of all 
stars with no Ha excess younger than 16 Myr and with masses 
in the range from 0.9 and 4M Q , marked as small dots in the 
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Left panel- 
from Siess et al. 12000) used to assign the age to the stars in the range of masses between 0.9 and 4 Mq 
are shown as thick dark dots and asterisks, alternatively. The other PMS stars are shown with large grey dots, whereas all other cluster stars are shown with 
smaller grey points. Right panel- Histogram of the number of PMS stars (i.e. stars showing Ho excess emission) per Myr as a function of age (solid line). The 
same histogram for all other stars is shown for comparison (dashed line). Note that for the these stars we have only sampled the age distribution up to 16 Myr, in 
order to avoid the significant field star contamination in the bluest part of the CMD. 



left-hand panel (note that the few remaining bona-fide PMS 
stars outside of these age and mass ranges, shown as thick 
grey points, and are not considered here). For the objects 
with no Ha excess emission we have only sampled the age 
distribution younger than 16 Myr, in order to avoid the signif- 
icant field star contamination in the bluest part of the CMD. 
For comparison purposes, their histogram has been normal- 
ized vertically (i.e. shifted by 1.1 dex) so as to match the 
distribution of the bona-fide PMS stars at the youngest age. 

The histograms of Figure|8] show that star formation in 
NGC 3603 has been ongoing for at least 10-20 Myr and no 
gaps are evident, at least at the level of resolution that we 
have adopted for the age (a factor of two, as shown by the 
size of the bins). If we consider stars with no Ha excess 
emission, for which no selection effects are present other than 
photometric completeness that is nonetheless always > 85 %, 
we would have to conclude that over the past ~ 16 Myr the 
star formation rate in this region has been progressively in- 
creasing. However, it is important to consider that many of 
the older stars might have migrated out of our FoV. Accord- 
ing to Rochau et al] (2010), the velocity dispersion of stars 
in the central regions of NGC 3603 is ~ 4.5 km s" 1 and ap- 
pears to be pretty constant in the mass range that they sample 
(~ 1.7— 9Mq). This implies that a 10 Myr old star would have 
had time to travel as far as 45 pc away from its birthplace, i.e. 
well beyond the ~ 5 x 5 pc 2 area covered by our observations. 
This might be one of the causes of the observed drop in the 
number of stars with increasing age shown in the histogram 
of Figure|8] and would also explain the somewhat different 
radial distributions of old and young PMS stars seen in Fig- 
ure [7] A survey of a wider area around NGC 3603 is needed 
to properly address the evolution of the star formation rate in 
this cluster. 

Interestingly, the age distribution of bona-fide PMS stars 



does not seem to differ in any systematic way within the er- 
ror bars from that of objects with no Ha excess, seemingly 
suggesting that the ratio of PMS stars with and without Ha 
excess emission is not a function of age. This might appear at 
odds with common wisdom suggesting that the efficiency of 
the accretion process at the origin of the Ha emission should 
decrease with time as a PMS objects approaches th e MS (see 
ISicilia-Aguilar at alj2010t iDe Marchi etaTll2010al and refer- 
ences therein). However, selection effects here can be impor- 
tant. In particular, as PMS objects approach the MS, both their 
Ha and bolometric luminosities decrease, but not necessarily 
in the same way. If the bolometric luminosity drops more 
rapidly than the Ha luminosity, PMS stars of older ages will 
become easier to identify for our method si nce it requires an 
excess emission at the 5 a level or higher (see lDe Marchi et ail 
l2010ah . 

This effect can have important implications on our under- 
standing of the accretion process and of the star formation 
rate, which we will address in a forthcoming paper compar- 
ing the star formation process in a number of young clusters. 
We have already started to study the stellar population in the 
30Doradus region in the LMC, applying the same observa- 
tional strategy described in this work, and using the recently 
released WFC3 observations of this region. By searching 
for objects with Ha excess emission, we have already found 
tantalizing evidence of multiple star formation episodes (De 
Marchi et al., in preparation), like in NGC 3603. Further- 
more, two distinct populations with ages of ~ 1 and ~ 15 Myr 
are present in the population of the star cluster NGC 346 in 
the Small Magellanic Cloud (De Marchi et al., in prepara- 
tion). Given the different environmental conditions and chem- 
ical compositions of the three clusters (Z Q for NGC 3603, 
~ 1/3 Z for 30 Dor and - 1 / 10 Z for NGC 346; see HA08, 
Andersen et al. 2009 and Hennekemper et al. 2008, respec- 
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tively), this similarity supports the hypothesis that continuing 
star formation could be the preferential channel for the forma- 
tion of stars in starburst clusters. 

According to Vinko et al. (20Q9|), the young cluster 
Sandage96 in NGC 2403 is known to positively exhibit a 
young population (10-16Myr) together with a relatively 
old one (32- lOOMyr) thus suggesting multiple star forma- 
tion events in a range of ages at least 4 times wider than in 
NGC 3603. A spread in the MS turn-off has been reported for 
clusters of intermediate age in the LMC and has been inter- 
preted as an age spread of ~ 300 Myr (see Milone et al. 2009). 
Moreover, the discovery of multiple stellar populations along 
the MS and red giant branch of a large numb er of Galactic 
globular clusters dPiottol 120081: iLee et al.1 120091) requires two 
or more bursts in the star formation history of the se objects, 
separated by at a least few 10 7 yr (see lCarretta et al . 2010, and 
references therein). 

Therefore, it appears that multiple generations of stars 
spread over a wide range of ages are present in star clusters. 
A detailed investigation will be necessary in order to under- 
stand what is at the origin of the observed age spread and, for 
example, what is the influence of the environment and chemi- 
cal/physical state (metallicity, turbulence, density, mass, etc.) 



of the parent molecular cloud on the formation of stars in clus- 
ters. 

Establishing whether age spreads like those seen in 
NGC 3603 are common in starburst clusters will have pro- 
found implications for theories of star cluster formation, for 
the meaning and determination of the IMF and finally for 
the general assumption that clusters are simple stellar pop- 
ulations. 
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